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SllfOPSIS 


The present work is based on two separately developed 
concepts of procurement and inventory operations. Por these 
concepts, two models are proposed. The first model deals 
with a multi-item inventory system involving single source 
and limited warehouse capacity. A dynamic programming 
formulation with usual cost terms is structured. Objective 
is to minimize the total cost per period. The total cost 
is expressed as a function of maximum inventory level. 

Depending on the situation under consideration, this variable 
can take discrete or continuous values. Therefore, both 
discrete and continuous dynamic programming have been used 
to solve the problem. 

The second model deals with joint replenishment of 
several items. The problem has '■ been formulated mathematically 
with the objective of minimizing the total cost. Three 
heuristic approaches have been developed. Heuristic one 
considers a planning horizon of finite duration. Heuristics 
two and three solve the same problem but the planning horizon 
is taken as infinite. For an n-item problem heuristic two 
generates n feasible solutions, and the feasible solution 
giving the minimum objective fmction value isselected. 
Heuristic three is an iterative procedure and it generates 
only one solution. 
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The three heuristics have been cempared for a set of 
45 problems generated randomly. It is found that heuristic 
three which is computationally very efficient yields better 
solutions in most of the cases. 



CHAPTER I 


IHTRODUCTIOH 


THe control and maintenance of inventories of physical 
goods is a problem common to all enterprises in any sector of 
a given economy. Inventories are maintained by manufacturing 
as well as non-manufacturing organizations. The economic 
significance of effective control of inventor 3 /' is without 
question very great, for the costs associated with procure- 
ment and holding of inventory are large in almost all Indus trials 
business and service organizations. It is no secret that 
inventories consume most of the investment capital in any 
enterprise. In fact inventories can represent more than 
50 percent of a company's invested capitaJL and as much as 
80 percent of its cost of goods sold [lo]» Because of the 
magnitude of costs involved and the widespread applicability 
of inventory problemss, it is obvious that inventory theory is 
a fertile field for cost reduction efforts. 

Al though literature on inventory can be traced back to 
the turn of the century, useful work has ensued only after 
1954. The later woiic attempts at structuring real life 
problems. The formulations, however, have been so complex 
and cumbersome that most of the research workers have concen- 
trated on only single product inventory models. 
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The present -work can he seen as consisting of two 
paits. The first part presents a model for developing the 
optimal procurement and inventory policies with respect to 
a cost function in a deterministic system composed of many 
items that may he procuied from a specified source. The 
system is subject to a restriction on total warehouse 
space. The policy variables for the system under considera- 
tion include the procurement level and economic order 
quantity for each item. These are determined with respect 
to the relevant costs, system parameters, and a warehouse 
restriction such that the minimum total cost results. The 
minimization is performed through the application of dynamic 
programming technique. The solution proceeds stagewise 
(where each stage corresponds to an item) with the aid of 
recurrence relations and a functional equation method. 
Computer programs for both dynamic programming (discrete) and 
dynamic programming (continuous) ha-ve been developed, sasi 

The second part of the presentation aims at finding 
the optimal inventory policies (i.e, when to order and how 
much to order) with respect to a cost function when it is 
possible to replenish many items jointly. In case of joint 
replenishment the economic lot size for each item cannot bo 
determined without considering the rest of the items within 
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the group. Bach time a hatch of product is procured,- there, 
is a common set-up cost for procurement in addition to the 
individual set-up cost for ea.ch of the items. Therefore, 
the economic order quantity for each item cannot he deter- 
mined without considering the rest of theitems which are 
jointHy rcp3.enishe d with the item in question. 

The economic incentive behind thedesire for joint 
replenishment as a group lies in the fact that if each item 
is replenished individually, it must he responsible for the 
full set-up cost for each of its procurement run. In a 
manufacturing context, economies may accrue throu^ joint 
replenishment as a consequence of savings in facility set-up 
and changeover time. In jointly ordering a number of items 
from a vendor cconomiosmay also accrue as a consequence of 
the preparation, processing and expediting of fewer oixSers 
and the receiving of fewer shipments. 

To solve such problems one common approach is to 
determine independently the lot size for each item as if no 
inter-relationship a.mong items exists, 8.nd then to coordinate 
the replenishment of the items by modif 3 ring the lot sizes and 
cycle times of the items. One approach used in practice is 
to decrease the cycle time (and, accordingly, the order 
quantity) of all the items with a low order fiequency, making 
it equal to the highest frequency item. The shortcoming of 
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tills approach, however, is that it is unreliable. The 
method does not always yield a solution with minimum total 
cost. In fact, it sometimes yields solutions which arc 
more costly than those in which all cycle times are simply 
left at their original values. 

In the following' chapters, three different algorithms 
are presented for determining minimum cost solutions to 
multi item problems of this type. The first algorithm 
assumes a planning horizon of a finite duration and the 
last two algorithms assume a planning horizon of infinite 
duration. Solutions are sought in which all demands are 
metE exactly and forwhich a reasonably good value of objec- 
tive function is obtained. 



OHIPTSR II 


LITBRATURB SIHW 


Inventory control is a field which is Kipidly growing 
in terms of available Iite3r8ture and an attempt to assess 
the worl which has been done already in this field and 
placing it in its proper 'perspective appears to be an uphill 
task. Keeping in view the objective of the present work, thi 
chapter concentrates on 'Multi-product Inventory Models'. 

An attempt has been made to give an ancount of the state of 
literature as it stands today. 

Sts.rr and Miller [23 ] appear to bo the first to 
introduce the concept of constrained multi-producb inventory 
systems. They assumed instantaneous rc’plonishmcnt and 
backorders wore not permitted. With the use of Lagrangian 
multipliois, an optimal sifcocking policy was determined for 
a given set of physical and economic constraints . Problem of 
this type has also been investigrted by Hansman [ 7 ], Buchan 
and Koenigsborg [ 3 ], and Hadley and I'Oaitin [ 6 ]. 

Haddor Ipy] considered a tvio product system and found 
two optimal policies. The first policy pertained to the 
joint ordering of products using the classical lot-sizo 
formula. In the second policy independent ordering of items 
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was peimitte d. Naddor and Saltzmann [18 ] considered a 
multiproduct system fjlaere tho objecti^^TC function was s sum 
of carrying and ordering costs. HoT-jevci"; the ordering cost 
was a function of items listed on a single order. The prob- 
lem was solved using numerical approxii'’ ations ana it gave 
subop timal results. Davis [ 4 J modified the TFaddor and 
Saltzman* s approach. He formulated an objective function 
that consisted of orderingj interest, obsolescence and 
holding costs. The holding cost tcrr.i included coots for 
warehouse space, labour and equipment. However, since the 
obsolescence, intci’est, and holding costs wore all of the 
same fo m (a percent of the va.lue of tho item) they were 
consolidated into a single term that wo.s similar to carrying 
costs considered in previous in'vestigCLtions . As before in 
the Haddor and SaHzman model, the terms wore independent and 
the total cost was simply the sum of costs for each item 
considere.d independently. 

Plossl [20] proposed a LIMIT system to find optimal 
ordering policies, once tho number of orders and reorder 
quantity are known. LIMIT is an acronym for lot -size Inventoiy 
Management by Interpolation Technique. It takes into account 
the practical limita^tion on the total number of sot-up hours 
which are avahlablo. Optimal order quantity is calculated by 
the ¥ilscn*s Lot Size formula and then the number of orders to 



satisfy tho demand are determined. 3y an iterative proce- 
dure the optimal mode ofordering is found for a given multi- 
product inventory p roll em. The method was designed to reduce 
tho total ordering costs. The method, ho^rover, did not 
recognise any product intera.ctions in tho form of limited 
storage spa.ee or iDudgetory restriction. 

Shu [22J has described a method for determining 
replenishment frequencies of two items jointly replenished, 
Sssentially, tho procedure developed a skipping rule for 
items. The skipped item is replenished ovciy K-th time (K ^1) 
that tho orders are placed. A set of graphs is provided to 
determine 2 and IT (whoro IT is tho number of times the product 
will be procured). If the number of items exceeds two, Shu 
suggested that the items be ranked according to decreaising 
order of sales and that the smallest demand item be chosen 
as skipped item. The same procedure is to be repeated until 
the value of E for the item in question ha,ppens to be 1 or 
there are no more items to evaluate except the largest demand 
item. Earlo [ 13 ] has modified Shu* s procedure. He su^osted 
that items be ranked in the decreaiSing order of domand/sot-up 
cost ratio and that tho item with smallest ratio be selected 
as skipped item. 
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Andrews and ilmiaons [ 1 ] have presented a branch and 
bound algorithm for finding the optimal ordering policy for 
a multi-product system with deterministic demand and no 
backlogging. Ordering costs were assumed to include a 
separate set-up cost *Kj^’ for each product ordered but 
witha a fixed savings, 2, each time all products were ordered 
simultaneously. The method, however, does not take care of 
the savings in cost when a few (not all) items are ordered 
at tho same point. 

Ignall [ll ] in his paper developed a method to 
determi-ne minimu:xL average cost ordering policies for conti- 
nuously revisited two product systems with joint set-up costs. 
Markov Renewal Rrog ramming was used to find the region in 
parameter space where a given policy was optimal. 

Herron and Hawley [ g ] have considered a multiproduct 
inventory problem where the objective function includes a 
term for storage anea cost which is a function of maximum 
expected inventory level G-crald F. Saxton [21 ] has considered 
extensively, prodiict interactions due to storage area. 
Saxton's work consists of developing an analytic method for 
deteimining an optimum cost policy as a function of order 
quantity for a total cost multi-product inventory model 
where demands are kno\m with certainity. The objective 
function is a sum of stora.ge, replenishment, carrying, 
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obsolescencG j area, bin, man-power and equipment costs. A 
computer simulation is included to varify input parameters as 
well as the resultant system cost as provided by the analytic 
model . 

lastly, Arvind R. lele [15] has proposed a model for 
a multiproduct inventory problem with, budgetaxy and storr, gc 
space limitations. Lele structures the situation as a non- 
linear programming problem. The problem is solved by lavidon- 
Rletcher-Powell Variable Metric Method with Golden Search 
technique applied for uni-diroctional minimisation. 

The review of pertinent literature in the field of 
inventory control has uncovered several aspects that merit 
further attention. One of these areas concerns the inter- 
actions which arc present in a miiLti product system. It x^as, 
therefore, felt that sdLutions making pa-ztial incision into 
this area would be an addition to this field. 

The first part of the present study dca.ls with a 
multi product inventorsT" problem where there is an upper 
limit imposed on total available storage space and items 
are competing for the floor space. The second part of the 
study investigates a problem where product interactions are 
present due to joint replenishment of two or more items. 



CHAPTER III 


TV70 PROBIEtlS; THEIR PORBRJIATIONS AI'TD SOLUTIOH 

PEDGEDHSES 

In this chCcptor we ccnsider two multi -item inventoay 
problems, and give their solution procedures. The first 
part of the chapter presents a dynamic programming algorithm 
to find out optimal ordering policies for a multi item 
inventory system subject to warehouse capacity constraint. 

In the second, joint replenishment of several items at a 
time is considered amd three heuristic jprocodures aie 
developed to solv3 this prctolem. 

3 . 1 Problem 1 ; 

A statement of the problem is given be lows 

5.1.1 Statement of the problem s 

A brief description of the decision situation under 
consideration is given below; 

Several items are mainta.ined in inventory to meet 
the demands. In satisfj/ing these demands, the systems 
manager finds it necessary to replenish li is stock for 
each item periodically. A limite d warehouse space is 
available for stocking the items. Tho manager must decide 
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wh.eii to procure oa,ch item smd how much of each item to 
procure so tha.t at no time the inventory level exceeds the 
maximum permissible by the available warehouse space. The 
decision environmentis composed of item dependent parameters. 
The objective is to dote mine the procurement level and the 
procurement quaiitity for each item so that the sum of aH the 
costs associated with the inventory system under consideration 
is minimized. 

3.1.2 Assumptions ; 

The following simplif 3 n-ng assumptions have been made 
in formulating the above problems 

1. The demand rate, replenishment mte and lead time for 
each item are known and constant. 

2. Backorders are p emitted. 

3. Item cost, procurement cost, holding cost a,nd back- 
order cost for each item are Iniovni and constant. 

4. Only a limited warehouse space is ava,ilable. 

5. The maximum inventory level in the warehouse is sum 
of the maximum inventory levels for individual items. 
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.1.3 Not e.t ions t 


The f oHovii g notations arc used to define the problems 


TSC 


TCj_ 

TG* 


HC. 

BG. 

PC. 


^i 


^i 

IGi 

^i 


^i 


Pi 


m 


i 



= total system cost per period 
= total cost of item *i’ per period 
= minimum total costof item i per period 
= holding cost of item i per period 
= backorder cost of item i per period 
= procurement cost of item i per period 
= unit cost of item i. 

= demand rate for item i 
= repleiiisliment rate for item i 
= totsi,! unit cost of d^ items 
= cost of holding one unit of item i in stock 
per period 

= cost of having one unit of item i bocklogged 
per period 

= cost of procurement of item i per procurement 
= numbor of periods per inventory cycle of’ ■ 
item i 

= load time for item i 
= Older quantity for item i 
= optimum order quantity for item i 
= procurement level for item i 
= optimum! procurement level for item i 
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¥ 


w- 


= maxiiaum inventoij- leTel in an in-ventory 
cycle of item i 

= maximum backorder level in an inventory cycle 
of item i 

= to tel a.vailable warehouse space 
= warehouse space consumed by one unit of 
item i 


5.1.4 Pr ob lem f o rmulation ; 

Since both demand and lead time for all the items are 
assumed deterministic , the resulting inventory process will 
be Qs indicated in Figure 1. The exhibited geometry of such 
a process will depend on procurement level, the procure- 

ment quantity, Qj_, the demand rate, d^, the procurement lead 
time, tj^j^ and the replenishment rate, r^^. Since these para- 
meters are item dependent, it is evident that the geometry 
will bo item dependent. 
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Math em a ti cal Re lat i onsh i ps 


Proin Figure 1, it is evident that the number of 
periods per invontoiy cycle for item i ares 


Also, the following relationships are ovidents 


■‘21^ <^1 


3i 4i'^ i 


^li "^ai 


Qi - + -^ 21 ^ *^^i " ^i^ 


^li ^ "^ai 


t^ . + t^. = 


I. + d.. tT, - L. 
1 1 Li 1 


From equations (a), (5) and (4), we get 


li = Qi (1 


) + Li - '^i-'fcpi 


The total time por cycle during which inventory is 
hold is (t2i + "t-^i). ILis is given by. 


*2± * '*^31 


h , h 

r. - d. 


From equations (5) and (6), we get 





*2i + "^31 = [Ql (1 - ^ ) + li - ^ 


¥e observe the f olloving rel ationship: 


(7) 




1 + 2- 

- di + di 


1 


d. 

d. (1 - — ) 
1 r.' 


( 8 ) 


From equations (7) and (8), we get, 

^21 '^31 [Qq (1 - ^ ) + • ■^ Li ] 


d, (1 - r^) 


(9) 


The total time for which a backorder exists is 
(tfi + ■fc^jL)* This is given by 2 

B. B. 

■^li ■^4i = ^ ^ 


( 10 ) 


But, since B^^^ = “^i * "'"Li ” ^i» 


+ t^^. - (dj_.tj^j_ - lj_) „ g__ + ) (11) 


li 41 


1 r . 

ill 


Prom equations (8) and (ll), we got the following 
re lati ons h ip ; 

"^li '^4i "" ^^i* ^Li ” ^i^ d7 

di (1 - 


( 12 ) 
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Total Systen Go, 31: ; 


The total cost per period for each item will be 
a summation of the iton cost per period, the procurement 
cost per period, the holding cost per period and the back- 
order cost per period, that is; 


TC^ = ICj_ + PCj_ + + BO^ (15) 

The item cost per period will be the product of tho 
item cost per unit and the demand rate in units per period, 


10. = C. d. (14) 

Tho procurement cost per period will be the procure- 
ment cost per procurement divided by the nunber of periods 
per inventoiy cycle, 


p. .d. 



(15) 


Holding cost per period will bo tho product of the 
holding cost per unit per period and the avc-rage nambor of 
units on hand during the period, 


HGj_ 


1 

-2’ 


n^ i'^2± ^ ^ 


li: 


m. 


Substituting for 
we get : 


(^21 + h 


and in the abo\e 


c quation, 



i 


HO^ = 




c . 


d , 


2Qi (1-5^) 


[Qi (1 - - «i.t^,.J 


(16) 


Backorder cost por period will Tdo tko product of tlic 
backorder cost per urit per period and the avera@5 number of 
units backordered durmg the periods 


®°i-= 2 




m. 

1 


0 r? 


BOi 


2Qi (1 - ■— ) 


(17) 


The total cost per period will be a summation of the 
four cost components given by equations (14), (15), (16) 
and (17). Therefore, 


TCj_ 




£ili 


h. 

± 


(1 - ^ ) 


[Ql (1 - ji) H- li - ^ 


2«1 ^ ) 


(18) 


The tota.l system cost will be given by. 

**. 

n 

TSG = TGj_ 


(19) 
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Warehouse Oapacity Gonstraln.'t ; 

The i-th item in the inventory system cons ■ones certain 

amount of x^arohouso space j 'There exists a certain amount 

of total warehouse capacity; k¥. The maximum accumulation 

of inventory for the i-th item, will consume w^ 

units of wa.rGhousG space. Therefore^ the restriction 
n 

^i^i ^ ^ must not he violated. In the follox'^ing 

i=l 

section, a dynamic programming algorithm is presented that 
will find the optimal procurement and inventory policies in 
the face of this restriction. 

3.1.5 Solution procedure ; 

The first stop in finding tho optimal procurement 
and inventory policies for the inventory system, is to 
determine whether tho constraint is active or not. Below, 
we give tho procedure for finding the optimal policies when 
tho constraint is active as well as when it is not antive. 

1. Optimal Policies xjhcn Constraint is Inantivo; 

The minimujo cost procurement level and procurement 
quantity for each item may he found hy setting the partial 
derive, ti'ves equal to zero and solving the resulting equa- 
tions. Prom equation (18) we obtains 
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P, . 


h-.Q. (1 




+ + 

2Q. (1 - -i ) 2Q. (1 - -3. ) 

X ^ ^1 

(20) 


Taking tke partial derivative ifith. respect to Q^, then 
with respect to “ ^i^ setting both equal to zero 


give 


^TC^ 


a . 

- h)" 

0 2 2 " o <3. 

2Qj_" U - ^) 
i 


t. (a..t^. - ip' 

P d. 

2«1 (1 - V > 


( 21 ) 


- h) 


h^ + 


(-'rhi - 

Qi {1 - ji ) 


Ij, (dl-tj,. - ip _ 

ar - “ 

Qi (1 - ) 


( 22 ) 


Equation (22) can be written as; 


. tx • — li- 

1 Li 1 


^i “ r7^ 
h. + b. 


( 25 ) 
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Substituting egustion (23) into (21)^ fre get 

2p, .d. 


qJ 


= [ 


■X '1 


d. 

(l - ) 

^i 




(24) 


Substituting equation (24) into equation (23) gives; 

d^ . 


* 

^i 


= d-.t^.. 


- [ 


• d. (l T* . 1 

i IT 

— ]‘ 

b, (1 + — ' 


(25) 




; 


Equation (24) and Equation (25) nay now be substituted 
back into equation (20) to give an expression for the ininimurn 
cost, if tor seveial sinplying stops, we obtains 


* 

TO. 

1 


= C. 


r‘*i+ C 


2 (1 - ^) Pi hi b^ , 

^ ] ( 26 ) 


(hf + h^) 


2, Optimal Policies when Constraint is Active s 


Since I. consumes warehouse snaco it is necessary to 
1 

express TG^ for each valueof over a reasonable range. 

These values make up cost functions for thedynamic programming 
algorithm. 

Equation (5) msy he solved for - l^^) 

giving; 




d. 

Qi (1 - “ ) 

i 


^i 


(27) 


Substituting equation (5) snd equation (27) into equation 
(18) gives: 
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1C. 


p. .d . h . . I . 

P] . li 1, . 1 1 


2Q. (1 ~ 

1 

V [Q, (l-^>) - I,]2 


2Qi (1 - ) 

1 


(28A) 


or 




p . • d . h. .1 . 

+ -Q— •*■ a-- 

2Q,- (1 ~ •— ) 

i 

d . 

tl (1 - ) 

+ ^ i Q. - b^.I. 


b, .1, 

1 1 


d. 

2Q, (1 - ji) 


( 28B) 


Taking the partial derivative of TG^^ with respect to 
Q.^- and sotting equal to zero we gets 




Pi-di 

^ 

^i 


hi V 

arp’ 

2Q.^ (1 - ^) 

X ^i 


d. 

bi (1 - ™ ) 

X X 


b. .1/ 

3 - 1 


2Qi^ (1 - ~") 

i 


0 
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This givGs; 


Q- 


1 


d, P 

2 p. ,d^. (1 - ^) + I ^ (b + bj 

i 


a . 

(1 « -2.) 


1 1 

— 


0 . 
j_ 


29) 


From equations (5) and ( 24 ) ^ we get 


2 Pi-q (1 - + b" {hi + bi) 

Lj_ ^ -j- - [ — ■ ■■ ] 


'■i* "Li 


b. 

1 


(30) 


Equations (29) and (30) give tho minimum cost and 
tlie minimum cost as a function of and other parameteiB» 
The minimum cost may be expressed as a function of and 
other parameters by substituting (29) and (30) into (28A) and 
modifying the last term. This giiresj 


TCj 


= C d + 4- 

^ 0* 

“1 


h. I.‘ 
1 1 


* d . 

2Q. (1 - ^ ) 






b. (d..t^. 
1 ^ 1 Li 


1* )2 


d. 

2Q* (1 - ^ ) 


(31) 


^'quations (29), (30) and (31) maj?- now be used to 
develop cost functions needed by the dynamic programming 
algorithm. In the dynamic programming algorithm presented 
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here, each stage corresponds to an iteni. Let g^ 
represent the total cost when inventoij'' level is and 
Ic^ fj_(¥) represent the minimuia cost value at the i-th stage 
(i,o. i-’th item) -vdaen waie house space available is *¥' 

(¥ is the state vci.riable) o 

The solution proceods stage-wise with the aid of 
I’ecurrence relations technique, First, the cost expected 
from the fixot stage (iteiu) if all warehouse space is 
a.llocated to it is determined frcm f-]^(¥) = g^(l 2 _w^). 

The compatoHons .for fp(¥) are completed and results are 
entered in the first stage of the solution table. In 
general, for the j--th stage, following recurrence relation 
is usedj 

f.(¥) = Min [g. + f. n(¥ - 

Od-w.-^W 

(32) 

The results for each stage are entered in the solution 
ui'blo. The solution table is iiicnused to find the optimal 
policies forchis constrained inventory system. The minimum 
total system cost is searched in the last sta-gc of the solutioiq 
table and corresponding inventory level is noted. Then we 
trace bach through all the stages and invuntory level at 
each stage (item) are noted. The optimal values of invent oiy 
levels so found i.o., Ij_’s, along witJi the equations (29) and 
(30) can be used to find out optimal order policies for this 
inventory system. 
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5*2 Problem 2 % 

The problem formulation and solution procedures for 
a multi-product in'ventoiy system with joint replenishment 
of piroducts are presented in this section. Three heuristic 
proceduies have been developed to handle both the cases 
when planning horizon is fini-fe as well as when it is 
infinite . 

3.2.1 Nomenclature ; 

The notations used in this section are given below^ 

^ = total number of items 

N = total nximber of set-ups in the planning horizon 

(for heuristic 1) 

= total number of set-ups por year (for heuristics 2 and 3) 
N* = economic number of set-ups per year 

= number of set-ups for i— th item in the planning 
horizon (for heuristic l) 

= number of sot -ups per year for i-th item (for heuristics 
2 and 3) 

N* = economic number of set-ups for item i in the plaming 
1 

horizon (for heuristic l) 

= economic number of set-ups per year for item i 
(for heuristics 2 and 3) 

= relative ordering freg.uency- for item i (kj_ = 
k* = economic value of k. 

T 1 
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rolati'TO ordering frequency for item j when item i is 
ordered at each set-up 
corrected value of k.. 

relative ordering frequency of item i at the r-th 

iteration 

corrected value of 

planning horizon 

number of periods in the planning horizon 
number of time unite per year 

tinie duration in an inventory cj^cle for which a, unit 

of item i is backlogged 

order quantity for item ’i' 

economic order quantity for item i 

annual demand for item i 

demand rate for item i 

fixed cost of talcing a sot-up 

variable cost of taking a set-up for item i 

cost of holding one unit of item i for unit time 

cost of having backlo^ed one unit of item i for 

unit t ine 

warehouse space aoccupied by one unit of item i 
total available warehouse space 
cost of occupying one unit of waiebcuse space 
during the planning horizon (for heiiristic 1) 
annual cost of occupying one unit of warehouse 
space (for heuristics 2 and 5) 



YC. 


TC 



variable cost for item i in the planning horizon 
(for heuristic 1) 

annual variable cost for item i (for heuristics 2 
and 3 ) 

toteil cost of the in'ventoiy system in the planning 
horizon (for heuristic l) 

annual total cost of the inventory system (for 
heuristics 2 and 3) 

annuc.l total cost of the inventory system when item 
i is ordered 8.t each set-up 


3.2.2 Case of Einite Planning Horizon s 
1. Problem Statement ; 

There aro n itoms stocked in a warehouse. The various 
costs, namely, 

(i) Fixed a,nd variable set-up costs, 

(ii) Holding cost, and 

(iii) Backorder cost 

arc known for oa.ch item. Decisions as to how many orders arc 
to be placed, when to place them and how much to order each 
time aro to bo made. Those decisions should minimize the 
total cost incurred over the finite planning horizon, T. 


2. Assumptions ; 


The model developed, is based on a sot of assumptions. 
The purpose of these assumptions is to simplify tho analysis 
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of the prohlem. Those assumptions aro descrihed below; 

1. There is a family of n products which can be ordered 
together. 

2. -An item is replenished at equal time intervals. 

3. Backorders are permitted. 

4. Rcplcaais hment is instantaneous for each item. 

5. There is only a limited stors,ge space available for 
stock holding. However, additional space can be 
rented o.t a constant rate. If there is some extra 
space it can be rented out at the same rate. 

6. The planning horizon is finite and known. 

7. Demand rate for each item is given end is constant 
ove r t ime . 

8. The demand for each item witliin the horizon period is a 
quantity greater than zero. 

9. There will be at least one set-up for each item within 
the horiz on period . 

10. The planning horiz oa consists of a finite n\imber of 
periods and a set-up can be taken up only at the 
commencement of a period . 

11. in order will be placed for each item at the beginning 
of the first period. 
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13. TiiG nuLabcr of in'ventctry cycles for each item, within 

the plaming horizon, is an integer. 

3. P rob lem Form ulat ion i 

The inventory process resulting from the problem as 
stated above is shown in Figure 2. 



Pig.2j Inventoiy Process for Problem 2 (Finite Horizon) 
Since the i-th item is replenished in equally 
spaced set-ups in the horizon period, the variable cost for 
this item consists ofs 

(i) variable set-up cost of set-ups? 

(ii) total stock holding cost in the plaming ho.il zon? 

(iii) total backorder cost in the planning horizon? and 

(iv) cost due to waie house space occupied. 


Thus, we can writs, 
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vo^ = i - -ti) 


+ i t?.dj^.iL .b^+ (Qj,-t^.d^) w^^.R 


(33) 


Tho total cost of the inventory system can bo written as. 


n n 

TO = ITS + ^ F. s . + 4- 7“ ( Q. -t . . d . ) X 

i^ 1 ^ 2 ^ ^W3_ ^ 3_; 


n 


- t.) N..h. + i ^t|.d^.F..b. 


n 

+ R (y^ (Q . -t . .d . ) w . - ¥) 
i=l 111 1 


(34) 


where the order quantity, is given by 

Td. 

Qi = y C35) 

To find the optimum value of t^^ which corresponds 
to minimum variable cost, substitute the value of in 
equation (33) and equate the partial derivative of YO^ 
wteth respect to t^ to zero . This gives , 

'w. .R + T.h. 

t. = -i — (36) 

%(h^+b^) 

Substituting this value of tj_ in equations (33) and (34), 
and replacing by Td^/y , we get, 
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VC^ 


■2N. 


1 X 3 . 


h.+lD, 
1 1 


+ R 


w^.d^ (2T,‘bj_ - w^,R) 

21T. (h. + b.) 

I X 


and 


n 


TG = 


T d. li..b. 

us + ^ 


i=l 


1 1 


w .d (2T.b - w..Pu) 

•*■ "" “X Th-Tb-T™"- J - 


(37) 


(38) 


4 . Solution Pro CO dure ; 

Heuristic Algorithm Is 

The solution procedure consists of tw parts. 

The first part consists of finding tlio values of for each 
item. In the second part, these values of aie plotted 
and all the points at which one or moro items are ordered 
arc counted. This gives the total number of set-ups, in 
tho planning horizen. 

The economic ordering frequonejr, for tho i-th 

item corresponds to tho minimum variable cost for the i-th 
item. The valxos of for i-th item can, therefore, be 
obtained by taking tho first derivative of equation (33) with 
respect to Hj_ and equating it to zero. This gives, 



31 


n. 


T - w^.R) 

2s ^ (li. + lj_) 


This value of is then corrected to ensure an 

integer number of invent oiy cycles in tbo plhinning horizon. 

A stop by step algorithm of the proposed hoiristic 

is given belowt 

Par t 1 s This part consists of obtaining the values of 

Step 1; Prom the givDn viluc of P find out the possiblo 
values which can take. Since the number of 
set-ups, for the i-th item should bo on integer, 

the possiblo values which can take will be the 
f ac to r5:i of P. I/ot these bo denoted by a^^, b- 2 * • “ ®'ps * • 
Thus, for each item i the number of set-ups for that 
item will be given by one of the a^'s. To illustrate 
this, let P = 12, then the values of a^ will be 
Op = 1, a 2 = 2, a^ = 3, a^ = 4, a^ = 6 and ag =12. 

So can take any of the values 1,2, 3, 4, 6 and 12. 

Step 2; Find out the va-lue of for i-th item frera 
equation (39)* 

Step 3; If this value of ]Ij_ is less than one then put = 1 
and go to Step 9, otheivrise continue. 

Stop 4! If the value of calculated in Step 2 is greater 
than P,thcn put N* = P (This is done because there 
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cannot be more than P set-ups in the planning 
horizon) and go to Step 9 othervri se continue. 

Step 5 5 If the -value of is equal to one of the factors 
of P, say Up 5 then put = a^ and go to Step 9j 
otherwise continue. 


Stop 6 5 


If the value of H. lies b e-tween b.. and a-,T thon 

1 J — 

! 

put 11^ = a^. and compute the value of -variable 
cost VC il’ i-th i-tom by the following 

formula; 


T .d- h..b. 
’^^il ^ ^i*®i hT+'b^ 


w..d . (2T.b . - ¥ .R) 

^ - , - - - 

(h^+b^) 


Step 7: 


Take the next higher value of a i.e. sind let 

^i ~ °'J-^-l’ compute the variable cost 

wi-fia this value of N?. let it bo V’C^ 2 » 




+ 


T^.dj_ hj.b^ 

2ir" h.+b. 
1 1 1 


+ R 


w^.dj, (2T.bj_ - Wj_.R) 

2< (hi + \) 


If the value of VC is less than or equal to 
VC ^2 \ = ^i SO -bo Stop 9, 


Step 8: 
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otherwise put = n". 

Step 9s Repeat Steps 2 through 8 for all itenB. If the 

values of for all the items have heen calculated 

then go to Step 10. 

]p This part of the algorithi? finds out the total 

number of set-ups, N, thctotal cost of the S 3 ^stem 
and the economic order quantity for each item. 


Stop 10: Using the abovo calculated values of find 

out R. for the i-th item X'rh ore R. = P/R.. 
u 1 ' 1 


Stop 11s An order for ovoiy item is placed at the first 
period and succossi-vc orders are placed after - 
every Rj^ periods. Rind out the periods at wleich 
an order is placed for tho i-th item, let these b o 


denoted by p^^x • • • :> Pj_p (where the 

maximum value of p^^ can be P and minimum val.uc 
can bo l). 


Stop 12: Repeat Stop 10 and Step 11 for c^-ll tho items and 
find out the values of Pj_i.*''5 for all the items. 

Stop 13 : Count those- periods at which one or more itoms 

are ordered. The total number of such periods give 
the tot ill number of set-ups, IT. 

Stop 14: Obtain tho economic order quantity for i-th item 
fron the following formula: 



stop 15 J Obtain the total costof tho system by substituting 
values of iT and ilf 5 in the equation ( 38 ) . 

3.2.3 Paso of Infinite Planning Horizon ; 

In this section two procedures are described 
which deal with a joint replenishment problem of a. multi - 
item in ■vent ory system when the planning horizon is infinite. 
The statement of the problon is gi-ven below? 

1 . Problem State mo nt ; 

In the previous model we treated the planning 
horizon as finite. The same problem is dealt here with an 
infinite planning horizon. The decisions are to be made 
such that the total cost per year is minimum, 

2 . Assumptions s 

The f iiBt five assumptions made in the previous 
model hold true for this case also. Following additional 
assumptions have been made; 

1 . The set-ups are taken at equal time intervals. 

2 . Planning horizon is cither infinite or it is 
so long compa.ro d to the C3rclo time that it 
can botaken as infinite . 



3. Atloast oric item is ordorod at every set-up. 

4* Total number of sct-up.s per year need not bo an 
integer . 

5. All items arc ordered at the first set-up. 

6. Annual demand for each item is a quantity greater 
than zero . 

3 . Problem Pormulation r 

The inventory process resulting from the problem 
statement given above is sho'wn in Figure 3. 



Fig. 3 j Inventory Process for Problem 2. (Infinite 
horizon). 

If the i-th item is replenished at every 
set-up, then there will bo set-ups por year for this 

item. Hero, k^ must bo a positive integer and is called as 
relative ordering frequency for the i-th item. The annual 
variable coot for the i-th item is given by; 
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YC^ 


w t . ,D. 

JL g. + 1 /n _ . 1 1 

kj_ '’i ^ ^ U.j_ 


mk. 

) (-# - 1-- ) 


IT 


1' k. 




2 t . , D . 

^ -S- '’i + (Qi--Vi) 


m 


(40) 


Ths ijctsl ciriinial cost: foptihB invon'ijor'y system is given by 

n -^T n 


TO 


NS + 


i=l ^i 


w ij. t . . D . mk . 

^ s.. + 4 . -^)(-^ ~ tj 


k, 


i=--l 


m 


li 




n 


2 ®1-® 


1 i£^ ' “ "i ^^E7 “i + “i («i- 


where the order quantity, Q. for item i is given by. 


'^i = 


D . . k . 

1 2 . 

N 


(41) 


(42) 


In the expression for YO^, substitute tho value of 
Qj_ from equation (42) and equate the partial derivative of 
YO^ with respect to t^^ to zero. This gi'ves optimum value 
of tj^ as; 


Pxkr-/ . + .h. 

t.. = —3—3 1 — 1 

^ N(h^ + bj_) 


(43) 


Substituting this value of t. in oqua,tions (40) and (41) 


ar3d replacing Q. by D-.k^/N^, we get; 


1 
Ns 


YCj^ 


D. .m 
_i + ..1— 
kj_ ^ 2N 


h. ,b. 
-JL-i ir 

hi+bi'^'-i 


w. .D. .k. (2m. b. - E.w. ) 

^ ^ 2m4 (h. - 


( 44 ) 
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and. ^ 


TG 


JL D u h..b. 

2 . C''’p“ + ~7Itt“ • ■u^V^ k + 

i=l ^ 


+ „ (aa.b, - BwJ, 

sti'n (hj_ + b:i •'- 


R,W 


(45) 


G 0 lu't ic n i 'TO c u dura t 


Huuristjic Ai^^oritim 2: 


In tfic ;ibovr'; problora, two policy decisions are 
to bo t-'jJcon: t,!iv tot'il nuabor of sot~upts per year F and 
the nur.ibci' of act -ups for individual itens F/k^. Dor a given 
value ol F, iho tcKi FS in equation (45 ) is constant. The 
valUw of k| { oi’ the i-tli itci can then bo found by differen- 
tiating oqu.' lion (44) partially with respect to k^^ and putting 
it wqual tu zujx) , This gives? 


N [ 


2®!-“ + '’!> 


^.1 •b* 4- lUw. .D. (2mb. - R.w, ) 

i. X 3L 1 1 


]” (46) 


or. 


F 

IT” 

1 


D. .h. .b- + R.w. .D. (2inb.-R.w. ) 


1 1 11 ^ 1 1' 

20j_.n (h^ + b^) 


V (47) 


V/e have assumed that atleost one item ‘will be ordered 
at 0£!ch set-up. This means at least one of the k^^'s will be 
equal to 1, The procodire adopted here can now be explained 
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as follows; First find out tho values of for each item 

from equation (47). Since one of the k^^'s should ho 1, 
therefore, start with, the first item and put tho correspond- 
ing valuLG of k^ equal to 1. This gives a particular value 
of N. From this valuo of I find out tho values of kj_ for all 
other items. Calculate tho total cost of the system with 
these values of F and hj_'s. Kext for the second item put 
k^ = 1 and repeat tho whole process. The above procedure is 
followed for all it3ms and tho one which gives minimum total 
cost is selected. 

A stop hy stop algorithm for the ahov3 procedure is 
given below; 

Stop 1: Use the expression (47) to find out the value of 

U/kj^ for tho i-th item. Let this value of U/k^ be 
called as Repeat this process to determine 

values of for all items. _ 

Stop 2: Atlcast one of the items is to be ordered at ctll tho 
setups or in othor words it can be said that atleast 
one of the k^'s should bo one. let it be for tho 
i-th item i.o. , k^^ = 1. So the value of N (i.e, total 
number of set-ups) will bo tho same as 

Step 3; If for any other item j ( j ^ i) the valuo of F . is 

J 

greater than or equal to the value of IT calculated in 
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Step 4 


Step 5 ; 


Step 6 


Step 2 then for that item put . =1 aiic. go to 

<3 ^ 

Step 9y otherwise continue . 

Galcula-te the value of k.. for those items only for 
which Uj is less than The rolue of k^.^^ for such 

an item can he calculated from the relationship 
given below; 

3i. 

k . . = 

11 N . 

J 

If for any item j, the value of k.. is an integer 
then put = k^.. and go vo Step 9, otherwise 

J -k tJ 

continue. 


If k^.j^ is a fraction then an integer value of it is 
to bo taken, which gives the minimuin vr.riablo cost f o r 
•the 3 -th item. In such a case k.. can talce either 

J 

lower integer value, soy or the newt higher 

integer value, say take the lower integer 

value and calculate the vo.ria.ble cost, soy 

t. — 

for the 3 -th itoia with the value of IT calculated in 
Stop 2. The •'/a.luo of con bo calcula.tod by 

using the following equation; 


U. . s . D . ,m 


YCj, 


JL -j- i 


■^1 


h ..b . 

_J 1 V 

2Hi • h^-.b. • ^3il 


R.w ,^ (2in.b.. - R.w,.) 


2ra.irp(h. + b.) 

J J 
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Step 7: 


kj.^ 2 i' again 


Take th.e next higlier integer value, 
calculate the variable ' cost, say 
same item, using the following equation? 


YG 


32 


^ji2 


I) . .m 




k... 


2li^“ * h^+b^. • ^oi2 


R.w-.D-.k-.p (2m b. - H.w.) 

-j. „„ J ti — 

2m. (h. + b.) 

1 J J 

Step 8? Gompaie the values of "70.^ and 70-2. It ^1 

c 2 

greater than or equal to 70^2 then put k^^ = k^^ 

c * 

and go to Step 9, othcr>a.sG,H put k^^ = 

Stop 9? Repeat Steps 3 through 8 for all tho values of 3 
(i.c., for tho value of H calculated in Step 2 
find out the values of k*^. for all tho items). 

c 

step 10; Oalculato tho total cost for the set of 

R values calculated abovo. The total cost, TC^, can 
be computed using the followiag equation? 


TG. 


n 

N.S + >" 

i=T 


c 


K.s. I)..m h..b. 

1 .1 . _J ...l.'.- J. . k° . 

c ^ 2U. • h>b7 31 


k7. 

31 


3 


R.w..D..k'5. (2mb. - R.w.) 

3 — 3 — ii 3 i-]_ R.w 

2m (h^ 


Step 11? Repeat tho Steps 3 through 10 for i - 1,2,..., n 
(i.o. putting kj_ = 1 for the next item calculate 
the values of R and k^j_’s and corro spending total cost). 



41 


Step 12; Select that value of 1 and set of k°.'s for which 

t3 ^ 

the total cost TC5j_ is minimejn. let the selected 
values of . ’s be represented as 1^. * s for all 
items . 


Step 15 i 


Oalculate the eccnomic order quantity for all the 
items using the following expression; 


Heuristic Piece dure 3; 

The method described in this section is an iterative 
procedure. To start with, a set of values of k^^ aie arbitrarily 
selected and the value of H is determined. The corresponding 
value of the objective function is then calculated. The value 
of the objccti've function is improves d at each iteration. If 
for any two consecutive iterations, value of the objective 
function remains the same, then the procedure is stopped. 

The logic behind the proposed procedure is- as follows; 

For each set of values of k^^'s there exists a local 
minimum. The value of H which comes pends to this minimum 
can be obtained by taking partial derivative of TC, in 
expression (45) wiiii respect to N and equating it to zero. 

This gives; 
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M 


2 . 
m 2l 


n 


D,- . 


h. ."b. 

1 1 


4- 

T 


i=l 


n D . . w . k . ( 2mlD . -R ,w ) 

k. + R7-i^=L 1 1 

i=l 


{ h . +b . ) 

1 X 


u. 


2m 


( S + 


n 


i=l 


fi- 


J 


(48) 

Again, a local minimum exists for ovqtj valao of IF an{3 corros- 
ponding values of k s for all items. For a given value of II, 
the term IIS in the expression for total cost is a constant. 

The optimal value of k^ for item i can, therefore, ho 
determined by talcing the partial derivative of variable cost, 
VC^, -wlih respect to k^ and equating it to zero. This gives s 


k.. 


IT [ 


2si.m (hj^ + bj^) 


D- .m'^ h..b. + R.W..D. (2mb. - R.-w. ) 

1 IX XXX X 


1 

IT 


Tho economic number of se-Wups, IT, and rala- 1 

tivc ordering frequency k^ for tho i-th item can be deter- 
mined with the help of egressions (48) and (49). The value 
of IF depends on a prior knowledge of k^^ for each item, which 
in turn depends on IF. in iterative procedure is, therefore, 
adopted for determining II and k^ for individual items. 

For the r-th iteration in the iterative procedure, 

wc . 


(a) use the k. values obtained in (r-l)-th 

JL y X ""X. 

iteration when IF was equal to 0 -K.d 
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(b) Obtaia and each item. If the value 

inte^r, then take the lower 
integer as value of k. and calculate the 

1 X 

valuG of variable cost for i-th item, let it be 

'''^*^i,l* K’sxt, take the higher integer as value 

of P and calculate the variable cost, let it 
# ^ 

lower integer is talcen as -the value of k. 

i,r ’ 

otherwise, the higher integer is taken as k . 

1 5 P 

The iterative procedure will be toiminated in r-th 
iteration if k^^^ = for all items. 

The following steps arc involved in the proposed 
iterative procedure; 

Stop 1; For the first iteration put r = 1 and select arbitraiy 
values of k^^ for all items. In practice, it is found 
that values of for i =1,2, ...n, provide a faster 

convergence. (If the aihitrarily selected values of 
k . ’ s happen to be nea.r the optimal values of k. *s 

IP 

then the ccnv3rgenco with these values of s will 

bo faster. ) 

Stop 2; OalculaiB the total number of set-ups IT from the 
formula giv^n below, lot it be represented as 
Thus, we gets 
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Step 3; 
Step 4 i 


Step 5 


Stop 6 

Step 7 


o n ii . . b . 

>“■ L --i ir 4. -Q <r- "i'"ir 

i=i ‘ £3. “ (b7"+Tr7 


D^. (2mb^ - R. w^) i 


1 

1 T 


n 


s . 

1 


2m (S + YJ ~ ) 
1 =1 ir 


(50) 


Set T = r+1 and i = 0, 


Set i = i-'-l. Use the value of calculated in Step 2 
(i.e. in the previous iteration) to find out the new 
value of Ihe value of can be calculated from 

the formula given below s 




2 bp 


k. 


N 


r-1 


I). .la'ih . *b .4- R.w - . D. ( 2rab.. -R.-w. ) 

1 J_ JL 1 JL ‘ J,. i. 


]"=■ (51) 


If the value of k^^^^ calcuLaled in Step 4 is less than 
or equal to 1, then put k?^ = 1 (This is done because 
no sot-up can be taker at intermediate points) and 
go to Stop 10. If the value of k^^ is greater than 1 
then go to Stop 6. 

: If k?^ is an integer number then put k?^ = k^^ and go 

to StepilO, otherwise, continue. 


Take the lower integer as value of k^^, 


let it b o 


represented as k^^. CalciiLate the variable cost for 
i-th item from the formula gi-ven below, ^et it be 
70 Then is given by: 
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va., 

1.1 


vl • , -1 




. m h . . b . 

_1 1 k" 

2\-l’ ■ ii’ 


(2mb^ - R.w^) 


2m. R 


r-1 


(h, + b, ) 


Step 8; Take the next higher integer as valne of 


t? 


it be represented as k. . 

ir 

cost with this value of k, 
70^2 is given by; 


Calculate the varie.ble 
Let it be 




VC 


r~l' 


i2 


D^.m 

+ 


ir 


r~l 


h..b. 
1 1 

h.+b. 
1 1 


k!’ 

ir 


R.w^.D^.k^^ (2mb^ - R.Wj_) 

2m.iT^_^ (h^ + bj^) 

Step 9: Compare the value of 70^^^ with ^ 

c ^ 

less than 

Step 10 5 othervTise , put k?^ = 

Step 10: If the value of i isless than total .number of 

items n, then go to Step 4, otherwise j continue. 


Step 11; If k? = k. ^ f o r i =1,2,..., n (i.e. if the two 

consecutive values of k^^ are same for all the items) 
then go to Stop 12, other&rise, put k^^ = k?^ for 


1 — 1 , 2 , ... 


n and go to Step 2. 



step 12: ae economic number of set-ups per year N is equal 

out the economic order quantity for 
each item, using the following formulas 

D . . k . 


Step 15 


Qj = 


IT 


* 


Determine the number of set-ups for each item from 
the following formulas 


H* = 


IT 


ir 


i = 1,2, 


Step 14s SubstitutG the vclues of IT and kj_ fo: 

into equation (45) and calculate the total cost of 
the system. 



CHAPTER IV 


MBTHODOIiOGT TBSTIFG AED SCOPE FOR HJTHRE WORE 

In this chapter results pertaining to the performance 
of the proposed models on a few test problems are presented. 
At the end of the chapter, possible avenues for further 
work are su^- ost e d . 

4.1 Performance of Pvnamic Programming Procedure : 

Nine problems have been solved using discrete dynamic 
programming computer package and for each problem optimum 
ordering policies have been determined. The problem si 2 © 
is varied 'from two items to four items. For each problem 
size a set of three problans have been solved. The average 
computation time for these problems on the IBM 7044 Computer 
system are given in Table 4.1, For tho sake of illustra- 
tion, the complete results of only one problem, for which 
the data is given in Table 4.2, are presented. For this 
problem, the total available warehouse space is taken as 
16 units. The results for the optimum ordering policies 
(i.o. , procurement quantity and procurement level for each 
item ) are gi van in Table 4.3. 

The same problem has also been solved using continuous 
dynamic programming procedure. The computer packa^ by 
Mize [14] was utilized for this purpose. The optimum 
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ordering policies for the problem, as obtained by the conti- 
nuous dynamic programming prDceduie are given in Table 4.4. 

It is found that the discrete dynamic programming 
procedure takes a computational time of 2.2 seconds and. 
gives an objective function valueof 284.689 units. The 
continuous dynamic programming results in a better objective 
function value of 278.736 units but it takes 17.605 seconds 
of computational time. Therefore, it is inferred’- that 
continuous dynamic programming approach yields better 
results, however, at a higher computational time. 

4.2 Performancos of Houristlc Algorithms. ; 

The performance of the thiee heuristics discussed in 
Chapter 3 have been cem pared by solving about 150 problems 
of various sizes. The problem size is~ varied from 6 items 
to 40 items. For each problem size a set of problaas have 
been genera. ted randomly. Computer packages have been developed 
in PORTRAY IV for the IBM 7044/1401 Computer system, for 
testing the proposed heuristics. 

The performances of the heuristics are evaluated based 
on cost index and average ccmputational time for a given 
problem sizo. The details of the two measures of evaluation 
arc gxvGn in follo’W’ing paragraplis. 
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Cost Index; 


For a given problem size, n, let TO and TO. be 

en in 

the total system cost using economic order quantity model 
and heuristic approach 'i* respectively. Then the cost 
index, is determined using the following relationship; 


Cl 


in 


TC. 

in 

en 


where, i := index of heuristic (i = 1,2,3) 

n = index of problem size ( n = 6,7,..., 40) 


For a given problem the heuristic which gives least value 
of the cost index is iaken as thc’bo^'-t heuristic. 


The three heuristics have been compared for a set of 
45 problems chosen randcmly. Heuristic one gives 'tho total 
system cost for a given planning horizon, whereas, heuristics 
two and three give annual total system cost. To compare 
heuristic one with heuristics two and throe, the planning 
horizon for heurislic one is taken as ono year. The results 
are presented in Table 4.5. The ro salts indicate that 
heuristic three yields least cost index valies in most of 
the cases (42 out of 45). In three cases the cost- index 
values obtained by heuristic two are least and in 11 cases 
these are very near to those obtained by heuristic three. 

The higher cost index values for heuristic one can be 
attributed to the following reasons which accrue frem the 
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tvro constraints invol-\ed in this heuristic; 

There is an upper limit on the total number of set-ups 
per year, 

2. The totel number of inventory cycles for each item 

has to be an integer. 

Further, examination of Tabic 4.5 indicates that the 
cost index values obtained for each of the heuristics is 
loss than 1. This implies that it i s economical to replenish 
several items jointly rather than ordering them separately. 

Average Gomputational Time ; 

For all the generated problems, the average caaputa- 
tional time is determined for each problem size. Average 
computational time values for different problem sizes are 
presented in Table 4.6. It can be seen from the table that 
for small size prcblems (upto 7 items) heuristic two takes 
tho least computational time. Heuristic three closely follows 
heuristic two. For large size problems (frem 8 items to 
40 items) heuristic three takes least time. A possible 
explanation for this behaviour of heuristics is given below; 

Heuristic two generates miiltiple solutions (equal to 
the number of items) for a given problem. Out of these 
solutions, the minimum objective function value solution is 
selected. Thus the computational time of heuristic two 
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incroasGs as th.G problom size incrcasos. Heuristic tharec, 
on the other .hand, is an iterative procedure, and as will 
he shown lattor, the number o f iteratiens do not depend on 
the problem size (the number of iterations, in general, are 
between 6 to 10). Therefore, for largo size problems 
heuristic three takes much less time compared to heuristic 
two . 

Based on the two measures of evaluation, it can be 
concluded that heuristic thoBe is the most efficient 
heuristic . 


Heuristic throe, being an iterative procedure, ic 
further tested to evaluate for its convergence. This is 
done on the basis of number of iterations it takes to solve 
a problem of given size and the percentage improvement in 
the total system cost. 


Bora problem of size *n', let TC^^ and TC^^^ represent 
the total system cost at the first iteration and at the last 
iteration respectively. Then the p ercentagems-improvement, 
PIj^, is given by the following relationships 


PI 


n 




,TC 


Ih 


X 100 


Table 4.7 gives the values of PI^ for 132 pi'oHlems 
ranging in size from 6 items to 4^ items. The second 
column of the table gives the number of pr|>l|l^m|,. g^l^^ 

eXHU -:: 

Acc, iuk-. 
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for the problem size under conside mtion. The third and 
fourth column of the table give the minimum andmazimum 
values of the number of iterations and percentage improve- 
ments respectively. As can be seen frcm the table, the 
number of iterations in most of the cases are less than 
10 (only in 4 cases the number of iterations exceeded 10). 

Thus, it can be concluded that number o f.i .'iterations do not 
depend on the problem size. Moreover, considerable improve- 
ment in the objective function value can be made in a few 
iterations. To illustrate it, for an 11 item problem 
percentage improvement of 88*13 is obtained in 9 iterations. 

In the previous paragraphs -we have compared the three 
heuristic s when the planning hbriz on is infinite. Now we 
shall evaluate the performance of heuristic one which considers 
finite planning horizon. The evaluation of this heuristic is 
done using the following procedure. 

The total system cost is calculated using the heuris- 
tic. The total system cost is also calculated when all the 
items are ordered at each period. The total system cost 
values using both the procedures are compared. 

Por a problem of size, n, let TG^ represent the 
total system cost when all the items are ordered at each 
period and let TO^^ represent the total system cost calcu- 
lated by heuristic one. The percentage cost saving, PS^, 
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obtained by usinga heuristic one is defined as; 

TC - TO, 

= ^ X 100 

An 

whcro n = index for problem size. 

Table 4.8 gi-ves the values of for 131 problems 
for the problem sizes ranging from 6 items to 40 items. 

Prom the table, it is clear that the use of heuristic one 
as cemparod to the proceduare of ordering each item at 
every period in the p laming horizon resulte in a saving. 

Por the selected set of problems, savings upto 97. 08 '-pe'rcent 
were obtained by the use of heuristic one. 

4.3 Scope for Puture Work ; 

In developing the methodologies presented in this 
dissertation, a few avenues for further work emerged. These 
are given below: 

1. In the dynamic programming procedure used it is 
assumed that the demand and lead time are deterministic. It 
will be worthwhile to extend it to the case where demand 
and/or lead time are stochastic. 

2 . In heuristics two and three it is assumed that there 
is no limit on the total number of set-ups per year. It 
would be interesting to analyse the problsn when there is 
an upper limit on the total number of set-ups in a year. 
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3. In the three heuristic procedures developed in this 

work, it is assumed that product inteiactions only due to 
joint replonishmont of the products are piresent. It would 
be worthwhile to consider the cases where additional product 
interactions due to warehouse capacity restriction and/or 
budgetary restriction are also present. 



Table 4.1s 


5S 


Average Computational Time for Problems 
Solved by Discrete Dynamic Programming 
Procedure 


S.No. 

Problem Size 

Average Computa- 
tional Time (Sec.) 

1 

2 

2.333 

2 

3 

3.075 

3 

4 

3.992 


Table 4.2: Data Used for a Tbree Item Inventory Problem 
Solved by Both Discrete and Continuous 
Dynamic Programming Procedures. 


Item 

No. 

Demand 

Rate 

Item 

Cost 

Holding 

Cost 

Backorder 

Cost 

Procu- 

rement 

Cost 

lead 

Time 

Space consxma- 
ed by one unit 

1 

6 

30.88 

0.30 

0.30 

18.30 

2 

4 

2 

4 

18.33 

0.24 

0.17 

17.50 

4 

2 

3 

1 

12,00 

0.12 

0.25 

15.50 

1 

1 


56 


Table 4.3 i Optimal Ordering Policies Obtained by 
Discrete Dynamic Programming Procedure 
for the Problem Given in Table 4*2 


Item 

No. 

Procure- 
ment level 

Procurement 

Quantity 

ffa^house space 
allotted 

1 

-14.08 

27.10 

4 

2 

-10.06 

29.07 

6 

3 

- 6.51 

13.31 

6 


Total Cost; 284.689 units 


Computer Time Taken* 2.2 seconds 


Table 4.4: Optimal Crdering Policies Obtained by 

Continuous %namic Programming -Procedure 
for the Problem given in Table 4.2. 


Item 

No. 

Procurement 

level 

Procurement 

Quantity 

Warehouse space 
allotted 

1 

-12.91 

27.26 

9.400 

2 

-11.20 

28.80 

3.200 

3 

- 8.80 

11.28 



Total Cost; 278,736 units 

Computer Time Taken; 17.605 seconds 


Table 4.5: Performance Comparison of Three Heuristics 
Based on Cost Index Values. 


SI. 

No. 

Pr oblem 
Size 


Cost 

Index Values 


Heuristic 

1 

Heuristic 2 

H-euristic 3 

1 

6 

0.482 


0.447 

0.331 

2 

7 

0.637 


0.610 

0.509 

3 

7 

0.904 


0.889 

0.887 

4 

8 

0.704 


0.636 

0.634 

5 

9 

0.723 


0.637 

0,636 

6 

10 

0.309 


0.288 

0.202 

7 

11 

0.482 


0.477 

0.420 

8 

12 

0.485 


0.478 

0.421 

9 

13 

0.421 


0.320 

0.249 

10 

14 

0,246 


0.230 

0. 216 

11 

15 

0.320 


0.294 

0.217 

12 

16 

0.591 


0.594 

0.592 

13 

16 

0,960 


0.V68 

0.944 ■ 

14 

17 

0.211 


0.187 

0,163 

15 

17 

0.318 


0.324 

0.283 

16 

18 

0.278 


0.249 

0,151 

17 

19 

0.591 


0.580 

0,555 

18 

19 

0.982 


0.953 

0.922 

19 

20 

0.527 


0.527 

0.486 

20 

21 

0.374 


0.363 

0.350 

21 

22 

0.723 


0.538 

0.487 

22 

23 

0.641 


0.623 

0,623 

23 

24 

0.623 


0.549 

0.556 

24 

25 

0.583 


0.548 

0.554 

25 

26 

0.712 


0.698 

0.488 

26 

27 

0.281 


0.299 

0.302 

27 

28 

0.136 


0.117 

0.115 

28 

29 

0.163 


0,146 

0.13T 

29 

30 

0.124 


0.114 

0.112 

30 

31 

0.143 


0.149 

0.129 

31 

32 

0.208 


0.187 

0.103 

32 

33 

0.189 


0.153 

0.153 

33 

34 

0.168 


0.131 

0.118 

34 

35 

0.178 


0.246 

0.135 

35 

35 

0.507 


0.532 

0.482 

36 

36 

0.126 


0.116 

0.112 

37 

36 

0.667 


0.675 

0.632 

38 

37 

0.292 


0.284 

0 , 169 

39 

37 

0.862 


0.831 

0.782 

40 

38 

0.224 


0,164 

0.103 

41 

38 

0.703 


0.675 

0.621 

42 

39 

0.152 


0.109 

0.099 

43 

44 

59 

40 

0.578 

0.173 


0.527 

0.162 

0.483 

0.160 

45 

40- 

0.80^ 


0.713 

0,629 
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Table 4.6s Performance Evaluation of the Three Heuristics 
Based on Average Computational Time. 


SI. 

No. 

Problem 

Size 

Average Computational Time (Sec.) 

Heuristic 1 

Heuristic 2 

Heuristic 3 

1 

6 

0.378 

0.267 

0.270 

2 

7 

0.420 

0.291 

0.295 

3 

8 

0.480 

0.368 

0.342 

4 

9 

0.503 

0.434 

0.320 

5 

10 

0.984 

0.483 

0.429 

6 

11 

0.729 

0.617 

0.392 

7 

12 

0.820 

0.723 

0.515 

8 

13 

0.793 

0.936 

0,581 

9 

14 

0.840 

0.927 

0.601 

10 

15 

1.082 

1.038 

0.642 

11 

16 

0.960 

1.026 

0.687 

12 

17 

1.020 

1.015 

0.730 

13 

18 

1.082 

1.061 

0.720 

14 

19 

1.146 

1.127 

0.761 

15 

20 

1.208 

1.201 

0.804 

16 

21 

1.403 

1.333 

0.785 

17 

22 

1.308 

1.293 

0.770 

18 

23 

1.382 

1.321 

0.920 

19 

24 

1.458 

1.431 

0.963 

20 

25 

1.489 

1.471 

0.916 

21 

26 

1.621 

1.599 

0.923 

22 

27 

1.672 

1.660 

0.978 

23 

28 

1.998 

2.016 

1.023 

24 

29 

2.123 

2.088 

1.042 

25 

30 

2.004 

1.950 

0.857 

26 

31 

2.087 

2.021 

0.885 

27 

32 

2.135 

2.080 

0.913 

28 

33 

2.263 

2.145 

0.947 

29 

34 

2.371 

2.212 

0.963 

30 

35 

3.398 

3.375 

1 . 021 

31 

36 

2.367 

2.341 

1 .028 


37 

2.428 

2.417 

1.056 

y c. 

33 

34 

38 

39 

2.491 

2.571 

2.483 

2.568 

1.073 

1.110 

y'~t 

35 

40 

2.648 

2.630 

1.138 
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Table 4.7: Convergence Evaluation of Heuristic 3, 


SI. Problem Ho. of Range of Total 

Ho. Size Problems Ho. of Iterations 

Solved 


Range of Percentage 
Improvements 


11 

r\ 

6 

3 

6 

— 

7 

28.34 


65,34 

d 

7 

3 

4 

- 

10 

22.84 

— 

38.85 

3 

8 

4 

4 

- 

8 

9.73 


59,83 

4 

9 

4 

7 


10 

55.27 

— 

83,97 

5 

10 

4 

6 

- 

8 

60.44 


64. 28 

6 

11 

4 

5 


9 

34.00 


88.13 

7 

12 

5 

6 

- 

8 

48.64 


58.13 

y 

13 

4 

6 

— 

9 

77.77 


80.23 

9 

14 

4 

6 


10 

64.66 


89.21 

10 

15 

4 

6 


12 

46.67 


65,15 

11 

16 

3 

7 

— 

9 

38.84 


57,36 

12 

17 

4 

5 


8 

43.86 


68.94 

13 

18 

5 

6 


9 

35.73 


58.75 

14 

19 

3 

8 

- 

9 

28,15 

— 

64,42 

15 

20 

4 

7 


9 

20.15 


52.34 

16 

21 

4 

9 


12 

52.54 


71,67 

17 

22 

3 

7 


10 

32.73 

— 

55.25 

18 

23 

4 

7 

— 

14 

58.63 

— 

84.58 

19 

24 

4 

6 

- 

8 

26.98 


63,72 

20 

25 

4 

6 

— 

14 

27.06 

— 

85.58 

21 

26 

3 

6 

— 

10 

52.93 

— 

88.90 

22 

27 

5 

6 


9 

54.79 

— 

88.63 

23 

28 

4 

6 


8 

29.36 


89.06 

24 

29 

3 

5 

- 

10 

26.63 


86.31 

25 

30 

4 

6 

- 

8 

29 .11 

- 

62,46 

26 

31 

4 

5 

— 

9 

26 .38 

— 

64.76 

27 

32 

4 

5 

— 

8 

9.22 

- 

35.27 

28 

33 

4 

7 

-- 

10 

42.82 

— 

64.17 

29 

34 

4 

5 

- 

8 

28.65 


68.48 

30 

35 

3 

6 

- 

9 

33.32 


47.46 

31 

36 

3 

9 


10 

37.40 

- 

53.58 

32 

37 

3 

8 

— 

8 

47.34 

- 

58.93 

33 

38 

4 

9 

- 

9 

50.05 

- 

69.28 

34 

39 

4 

7 


9 

51.04 

- 

73.82 

35 

40 

4 

6 

- 

8 

42.92 

— 

81,34 


SQ 


Table 4.8; Range of Percentage Cost Sa^vings Obtained by 
Heuristic One as Compared to the Procedure 
in which All the Items are Ordered at Each 
Set-up. 

(Total Problem = 131) 


SI, Problem Humber of Range of Percentage 

No. Size Problems Cost Savings 

Solved 


1 

6 

5 

52.79 


94.84 

2 

7 

5 

9.88 

— 

85.52 

3 

8 

4 

22.31 

— 

85.27 

4 

9 

4 

57.32 

« 

75.53 

5 

10 

5 

28.72 


94 . 76 

6 

11 

4 

11.72 

— 

72.87 

7 

12 

4 

19.94 

- 

76.93 

8 

13 

5 

38.32 

— 

84.50 

9 

14 

5 

16.99 

- 

87.73 

10 

15 

3 

58.20 

- 

82.47 

11 

16 

4 

67.37 

- 

94.08 

12 

17 

4 

68.62 

- 

85.68 

13 

18 

4 

32.64 


87.50 

14 

19 

5 

77.70 

- 

92.08 

15 

20 

3 

82.33 

- 

85.94 

16 

21 

4 

81.27 


82.82 

17 

22 

4 

48.82 

- 

79.03 

18 

23 

3 

32.74 


44 . 25 

19 

24 

3 

28.35 

- 

60 . 26 

20 

25 

5 

40.82 


62.35 

21 

26 

3 

70,76 

— 

97.08 

22 

27 

4 

44.26 

- 

76 . 56 

23 

28 

4 

21.79 

- 

53 . 74 

24 

29 

4 

52.35 

- 

74.86 

25 

30 

3 

21.59 

- 

42.73 

26 

31 

3 

20.95 


68.74 

27 

32 

3 

4.05 

- 

29.87 

28 

33 

3 

43.68 

- 

68.39 

29 

34 

3 

22.15 

- 

24.12 

30 

35 

3 

29.22 

- 

52.19 

31 

36 

4 

18,04 

- 

19.84 

32 

37 

3 

31.48 

- 

33.61 

33 

38 

4 

28.03 

- 

31.17 

34 

39 

3 

34.06 

- 

37.90 

35 

40 

3 

31.96 

- 

64.50 
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